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INDUCED GELATION IN A TWO-SITE SPATIAL 
COAGULATION MODEL^ 

By Rainer Siegmund-Schultze and Wolfgang Wagner 

Technische Universitdt Berlin and Weierstrass Institute for 
Applied Analysis and Stochastics 

A two-site spatial coagulation model is considered. Particles of 
masses m and n at the same site form a new particle of mass m + n at 
rate mn. Independently, particles jump to the other site at a constant 
rate. The limit (for increasing particle numbers) of this model is ex- 
pected to be nondeterministic after the gelation time, namely, one or 
two giant particles randomly jump between the two sites. Moreover, 
a new effect of induced gelation is observed — the gelation happening 
at the site with the larger initial number of monomers immediately 
induces gelation at the other site. Induced gelation is shown to be 
of logarithmic order. The limiting behavior of the model is derived 
rigorously up to the gelation time, while the expected post-gelation 
behavior is illustrated by a numerical simulation. 



1. Introduction. We consider a minimal spatial model of coagulation 
with the gelling kernel mn describing the coagulation rate and an inde- 
pendent motion of particles between two sites and 1 which happens at a 
constant rate for each particle irrespective of its mass. So the basic model 
is a continuous-time Markov jump process 

{Ct,m,j}, t>0, m = l,2,..., j=0,l. 

A state of the process is given by two sequences of nonnegative integers, 
both with only a finite number of nonzero entries, describing the numbers of 
particles of different masses at the two sites. As for the initial distribution, 
we consider the deterministic state given by the two length- 1 sequences 
{A^^oji {-^i} meaning a start with Nq,Ni monomers at the two sites. Denote 
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the total number of monomers by 

N = No + Ni. 

Given a state (Cmj), the jump rate (depending on the scahng parameter N) 
has the form 

iN{N-l) + KNY,U,j- 
The transition is defined as follows: Either, with probability 

m,j \ m,j 

a randomly selected particle (an aggregation of monomers) jumps to the 
opposite site, or with the remaining probability two monomers are randomly 
chosen and the corresponding particles coalesce provided they are not the 
same particle and they are at the same site. This means that an unordered 
pair of different particles of masses m and n at the same site coalesces at rate 
mn. Hence our model is a (minimal) spatial generalization of the Marcus- 
Lushnikov stochastic coalescent ML^^)(t) with monomeric initial condition 
[12, 13]. 

It is clear that we can describe the evolution as well as a colored graph 
evolution (with two colors): Start with Nq isolated red nodes and A^i isolated 
green nodes. Add at constant rate 1 an edge for each unordered pair of nodes, 
where edges between differently colored nodes are instantly removed as well 
as multiple edges, and let (maximal) connected components change their 
color at constant rate kN. 

We investigate the limiting behavior of this model for tending to infinity 
while 

(1) Ni/No ^ A, 

where A is a given parameter. The expected behavior is notably more com- 
plex compared to the well-known case where all particles live at the same 
site {mean-field model or homogeneous model). The mean-field model has a 
deterministic limit, described by a modified Smoluchowski system of differ- 
ential equations (with the modification being urged by the effect of gelation, 
meaning that, after rescaling time by t ^ t/N, at gelation time tgd = 1 a 
giant particle of mass comparable to the total mass suddenly appears in a 
phase transition and modifies the evolution). The more general model con- 
sidered here is expected to have a stochastic limit after gelation time, caused 
by the fact that the motion of (at most two) giant particles is no longer gov- 
erned by the law of large numbers. Moreover, a new effect of induced gelation 
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is observed: The gelation happening at the site with the larger initial con- 
centration of monomers immediately induces a gelation at the other site. We 
give a rigorous derivation of the limiting behavior up to the gelation time 
which is based mainly on the cutting techniques developed in [14, 15]. 

We emphasize that a mathematically rigorous derivation of the limiting 
behavior after gelation in our model can be expected to be quite difficult 
since the existing techniques either are based on a martingale approach and 
rely heavily upon the finiteness of the expected particle size per monomer 
[14], which breaks down at the gelation time, or they make use of random 
graph techniques [2] which admit no (or at least no immediate) generaliza- 
tion to the spatial model. Also, as will be seen later, the speed of convergence 
to the limit process is very slow, namely of logarithmic order in the total 
particle number, because induced gelation is an effect of logarithmic effec- 
tiveness. Hence, as far as post-gelation behavior is concerned (Section 4), this 
account is of mainly phenomenological nature and will deliver evidence for 
the expected phenomena by numerical simulation. In the rigorous part we 
use the notions "gelation time" and "gelation" in a weaker sense, meaning 
the explosion of the second moment. 

2. Limiting behavior of the spatial model before gelation. Here we study 
the two-site model at pre-gelation times. The parameter A indicates the 
dependence on the initial condition (1). 

Theorem 1. Let k> 0. Then there is a nonrandom time ig^j > 0, such 

that the rescaled quantities -/V~"'^Ct/Af,m,j our model with asymptotic pa- 
rameter A do uniformly converge in probability to deterministic quantities 
^tmj fulfilling a two-site analogue of SmoluchowskV s coagulation equation: 




,m—n,j 



-t,n,j 



.A 



0<n<m 



n>0 





.A 

't,m,l—j 



i + Kc: 



with initial condition 




A 



^0,1,1 — 



1 + A 



'-'0,m ,j 



,=0 



m > 1. 




m>0 

reach infinity for the first time. 



By a strong maximal positive solution of the system of ODE (2) we un- 
derstand a solution defined in an interval [0,t) which obeys the properties: 
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(a) < c^^ j for all ?n G N, j = 0, 1, t G [0, t ) , 

(b) the second moments (3) exist for t G [0,t) and fulfill 

rT 

J a^j ds <oo, T <i, j = 0, 1, 

(c) < t < +00 is the largest value with properties (a), (b). 

Theorem 2. The two-site Smoluchowski equation (2) has a unique strong 
maximal positive solution with i = t^^^. 

Next, let us consider the differential equation describing (as will be seen 
later, at the end of the proof of Theorem 2) the behavior of the second 
moments at both sites at pre-gelation times 

with initial condition 

Equation (4) is autonomous and extends in an obvious way the well-known 
equation a = valid in the pre-gelation regime for the mean-field model. 
Writing x{t) for (t^q and y{t) for a^i, we get the following system of two 
ODEs: 

(5) x = x^ + K{y — x), y = y'^ + K{x — y) 
with 

(6) xo = (l + A)-i>0, yo = A(l + A)-i>0. 

Theorem 3. /f k > 0, then there is a positive tgd < 2 such that x, y are 
smooth in [0,tgci) and we have x,y ^ +oo as t^tg^i simultaneously. 

Theorem 4. If k> and xq > yo, then x > y for all t G [0,tgei) and 
^ ~ (*goi ~ as t^ tgei while y ~ — Klog(tgei — t). 

Hence (except for k = 0) both second moments (expected particle masses 
at site j per monomer) (^tj^j = 0) 1) become infinite simultaneously (this is 
trivial for A = 1), an effect we call induced gelation in view of the mechanism 
lying behind this: If A 7^ 1, the immigration shortly before gelation time 
of large particles from the site with the larger initial number of particles 
(i.e., site in the case < A < 1) forces a gelation on the other site. The 
numerical studies with particle numbers up to 10^ show that this induced 
gelation exhibits some special features which distinguish it clearly from the 
"genuine" gelation taking place at the first site. In fact, for A < 1 the growth 
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of cj^o ^i^^ ~ ^gci) ^ while the growth of a^i is only of logarithmic 
order Klog(t — tgci)~^- Moreover, the proof of Theorem 4 indicates that, 
approaching tg^i, the growth of a^^ is finally completely caused by particles 
with large masses immigrating from site 0, while the "genuine" growth by 
coagulation at site 1 can be asymptotically neglected. 

On the other hand, it is clear from the logarithmic character of the induced 
gelation that it can be expected to be practically perceptible only for particle 
numbers even the logarithm of which is large enough. Take, for instance, 10^ 
particles at site and particles at site 1 with k = 0.3. As is known from the 
one-site model, at gelation time the size of the giant particle can be expected 
to be about 10(^/3)9 

= 10 . So in fact the mass of immigrating particles up to 
gelation time is bounded by this quantity. Therefore the maximum increase 
of the (empirical) value for o"^^ is most probably less than 6KloglO~4. 
That is an estimate of the maximal increase induced by site shortly before 
gelation happens there. The gelation at site 1 has to be accomplished now by 
coagulation of particles at site 1, since the immigration of massive particles 
from site rather rapidly stops after gelation, because these are rapidly 
consumed by the giant particle. The estimated remaining time it takes for 
site 1 to reach gelation is now about 1/4, meaning that in fact even for 10^ 
particles a perceptible delay occurs until the induced gelation actually takes 
place. Making k larger does not help too much, because then the tendency 
to reach equilibrium between both sites is stronger, such that it becomes 
harder to distinguish the induced from the "genuine" gelation, which would 
have taken place anyway at site 1, if migration was stopped shortly before 
tgci. Nonetheless, as we will see (interpreting the numerical studies), it is for 
the given choice of parameters not difficult to recognize the induced gelation 
by some characteristic properties appearing after tgei- 

3. Properties of the mean-field model. Consider the case k = 0,N = Nq, 
that is, the mean-field model. Then the graph is completely colored red 
for all times and we are in the situation considered in great detail in the 
monograph [2]. There the random graph model Q{N,p) is studied which 
describes a graph with N nodes, where for each pair of nodes an edge is 
present with probability p, completely independent from the presence or 
nonpresence of the other edges. Since the process described above connects 
an unordered pair of different monomers at constant rate 1 , the probability 
of the corresponding edge to be present at time t is 1 — e~*. Hence the 
distribution of {Ct,m,o}m>o is given by the joint distribution of maximally 
connected components of sizes m = 1, 2, 3 . . . in the random graph G{N, 1 — 

As is well known in both the coagulation interpretation (cf. the com- 
prehensive paper [1]) and the random graph setup [2, 5], the time-rescaled 
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(t t/N) model undergoes at t = 1 (i.e., at time in the old scale) a 
phase transition called gelation, which is characterized by the occurrence 
of a single "giant particle" of mass (,t^^^ = Q{N) [comparable to N in the 
sense that {^f^^"^)'^ = 0{N^^)] for t > 1. Up to this gelation time tgei = 1 the 
rescaled particle numbers N~^^t,m,o, m = 1,2, . . . , tend in probability to de- 
terministic quantities ct^m fulfilling the Smoluchowski system of differential 
equations 

(7) Ct,m = ^ '^('^ ~ '^)ct,nCt,m-n " mCt^m ""^i," 

0<n<m n>0 

which conserves the total (rescaled) mass 1, that is, J2m>o^^t,m = 1> t <1. 
At later times the limiting behavior is no longer described by this set of 
equations. As expected, the Smoluchowski system no longer conserves the 
mass for t > tgd, but it does not take into account the strong effect of the 
giant particle on the rest of the system. The actual limiting quantities of 
the stochastic coalescent differ from the solutions of (7), though sharing the 
property that total mass gets lost. 
The proper system of equations is 

(8) Ct,m = k '^("^ ~ n)ct,nCt,m-n " mCt^m- 

0<n<m 

Observe that in contrast to the infinite system of equations (7) the system 

(8) has the nice property that the finite subsystems obtained for m< M are 
autonomous. It is obtained from (7) by replacing the quantity J2n>o''^'^t,'n 
by 1, which simply means that particles are aggregated to others (including 
the giant one) at constant rate 1, representing the invariant rescaled total 
mass. 

The solution to (8) can be given in explicit terms for the case co,i = 1, 
that is, the monodisperse initial condition: 

(9) ct,m = t"'-^m"'-^e-'"'/ml. 

The convergence of the rescaled stochastic coalescent to (9), that is, the 
solution of (8), was shown rigorously [2] for all t>0. 

The quantities {mct^m}m>o, together with a mass of Ct,oo '■= ^ — J2m>o ''^(^t,r, 
assigned to the value +oo, form the Borel distribution with parameter t, that 
is, the distribution of the total progeny of a Galton-Watson process with 
Poisson offspring of parameter t. Observe that for t> 1 the branching pro- 
cess is supercritical and hence the total offspring is infinite with positive 
probability, leading to a mass defect: 

Ct,* '■= Yi ^^t,m < 1, 

m>0 

reflecting gelation. 
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Fig. 1. Plot of (t.t and at- 



While for the system (7) the second moment at := X)m>o "^-^'^i.m reaches 
infinity exactly at tgei and remains infinite, for the system (8) this second 
moment is infinite only at tgci (see Figure 1). In fact, at is just a reduced 
second moment, with the gel mass not being included. 

The second moment has a more immediate interpretation as a first mo- 
ment: It is the average mass of the particle where a randomly chosen monomer 
is aggregated into. So at gelation time this average mass becomes infinite, 
and it returns to finite values only because monomers "get lost" in the giant 
particle. 

Let us look at the behavior of the system (8) in some more detail. For this 
sake it is convenient to change from the consideration of particle numbers 
to the mass flow equation which is obtained if we consider the quantities 
Ct,m ■= nT'Ct,m. representing the mass concentrated in particles of size m at 
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time t. Now the evolution equation becomes 

(10) Ct,m = h"^ Ct,nCt,m-n " "T-G,m- 

0<n<m 

We introduce the ( probability) generating function 

u{x,t) ■.= 'Y Ct^^x"", 0<x< 1, 0<t <+oo. 

m>0 

Observe that due to the gelation effect at times t>l the total mass Q,* = 
u{l,t) is less than 1, hence u is for i > 1 no longer a probability generat- 
ing function in the strict sense. The evolution of u comprises all necessary 
information about the limit system including the mass of the gel-particle 
0,00 = 1 — u{^,t) and the expected particle size per monomer at = -^u{l,t). 
From a well-known formula for the probability generating function of the 
total progeny of a Galton- Watson process (cf. [8]) we derive the following 
semi-explicit formula for u: 

(11) x = ne*(^-"). 

Now observe some interesting consequences of this solution. First observe 
that for each t (11) has the trivial solution u = 1 for x = 1, and it is not 
difficult to check that for t < 1 this is the only solution of 1 = ne*^^""") 
fulfilling < n < 1, while for t>l there is a second solution with <u <1 
representing the true value of u{l,t) = Ct,* = 1 — Ct,oo', see Figures 2 and 3. 

Observe that the singularity at t = 1 appears as a cusp catastrophe in the 
full solution of (11), and the existence of the gel-particle seems to find some 
reflection in the solution u{x,t) = 1. 

So for t > 1 the total mass is given by the inverse of the function 

-logu 

It — > I = . 

1-u 

Also, it is easy to derive from this that the initial slope of 0,* at t = 1 is —2, 
respectively the initial growth speed of the gel-particle is 2, a fact mentioned 
in [2], page 123 (where it has to be taken into account that the time scaling 
there differs from ours by a factor of 2, such that his value is 4). Then, in 
course of time (t,* tends to zero exponentially fast. 

Another implication of (11) is a relation for the expected particle size per 
monomer at which reads as at = i^t(^ ■ From this it is easy to derive the 
asymptotics of cj in a right neighborhood of 1 as ~ (t — 1)~^ which shows 
that |t — 1|~^ is the asymptotics for a near t = 1 (for t < 1 this formula is 
well known to be the exact expression for a). 
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3 

Fig. 2. Plotofu{x,t),0<t<3,0<x<l. 

4. Expected behavior of the spatial model after induced gelation. The 

following remarks are only heuristical, but they are supported by the results 
of numerical simulations. 

It is clear that after gelation the spatial model will behave in a more 
complex way, since the limiting model cannot be expected to be nonran- 
dom. In fact, there will be two gel-particles at the two sites immediately 
after gelation, and these are not subject to a law of large numbers in their 
migration behavior. Consequently, the limiting behavior after tg^i must be 
expected to be characterized by jumps of the giant particle(s) from one site 
to the other. They immediately coalesce if they happen to be at the same 
site. The further time evolution of the quantities c^ mj strongly depends on 
the presence or absence of a giant particle at j. 

Hence, in contrast to the (mean-field) stochastic coalescent the (asymp- 
totic) gelation effect does not only imply a loss of total mass (caught by the 
giant particle) but additionally the spatial model shows up a stochastic lim- 
iting behavior. So the limit model cannot be described by a (deterministic) 
kinetic equation similar to (8). 

In order to describe the limiting object as a Markov process it is conve- 
nient to include the two gel-particles with their masses. This is not necessary 
in the mean- field model since the gel-particle simply represents the lost mass, 
whereas in our situation the distribution of the lost mass between the giant 
particles and their random locations are important. 




In order to unify the setup, we consider from this moment rescaled masses 
N~^m^t,m,j instead of rescaled particle numbers N~^^t,m,j- With this change 
of the point of view the Smoluchowski system valid for t < tgei gets modified 
to a mass flow equation (writing Ct',m,j '^^i,m j ) 

0<n<m n>0 

for < t < tgei- Let us denote the masses of the gel-particles by Ct)ooj,i € 
{0, 1}. They are zero for t < tgd. Finally, let TTj{t), j S {0, 1}, t>0 be two 
completely independent renewal processes with exponential waiting time 
(parameter k~^), with TTj{t) denoting the time since the last renewal event 
occurred. These two processes are the only random source of the limiting 
model. We get the stochastic evolution equations 

^ 1 aA aA 

0<n<m 

(12) 
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''^Ct,m,j I Ct,oo,j + ^ ^t,n,j I 
\ n>0 / 



'^^t,m,j + ^Ct,m,l-j 



and 



(13) 



+ l{7ri_j(t)<7r,(t)} 



A 



t-7ri_j(t),oo,j + Ct-7ri_j(t),oo,l-j 



+ 



for t > 0. This means in more rigorous terms that the hmiting process is 
a Markov jump process with piecewise continuous trajectories, where the 
jumps are determined by the (external) renewal processes 7rj(t), and between 
jumps the trajectory is described by an ODE (i.e., deterministic), while at 
the discrete jump times ti,t2, ■ ■ ■ given by 7rj{t) a gel-particle of mass Cti,oo,j 
changes from j to 1 — j and forms at site 1— j together with the gel-particle 
there (if there is one) a larger gel-particle of mass Ct^^ooj + Ct^,oo,i-j- Observe 
that we need to write the deterministic evolution equation between jumps 
as an integral equation since another gelation may have time to occur at the 
given site as long as Ct,oo,j ^^'^o with the consequence that discontinuities 
in the derivative occur. The trajectories of Ctmj piecewise differentiable 
between the jumps. 

So, from gelation time up to the first jump time happening after gelation 
[the pre-gelation jumps given by iTj{t) are irrelevant] the whole evolution is 
completely deterministic and may be described by (12) plus the equation for 
the gel-particle's masses (13) which can be rewritten as Ctooj = Ctooj^tj- 

Just as in the mean-field model after gelation equation (4) has to be 
modified, taking into account the strong effect of the gel-particle, and ceases 
to be autonomous: 



■=t,OD,j 

n>0 



It involves the second moment of the mass distribution, respectively the 
"third moment" in the traditional notation (referring to c instead of C,). So, 
while aQoo gives the speed of growth of Coo (of the gel-particle), the speed of 
descent of a now includes pQoo , where p is the second moment of the mass 
distribution, this being valid in the mean-field model as well as in our two- 
site model. The above differential equation for a may be expected to be valid 
except for the jump times (where right and left derivative are different) and 
except for the gelation time(s) (in our case there is more than one of them 
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with positive probability as explained above). In the latter case it becomes 
meaningless since a exhibits a singularity (and the indefinite expression • cxo 
is involved). 

In view of the fact that the consideration of the (probability) generating 
function of the mass distribution in the one-site (homogeneous) situation 
gives a detailed picture of the behavior of the limit model also after gelation, 
it is tempting to try this method also in our case to retrieve information 
about the deterministic behavior until the first jump happens. So one might 
hope for some more information about the asymptotics of Coo and a after 
tgci- Of course, it should clearly be expected that at the inducing site the 
asymptotics of a remains unchanged as in the one-site case \t — tgoi|~^ due 
to the fact that the influence of the induced site is negligible near tgei. It is 
less clear what will happen at the induced site in a right neighborhood of 
tgei, since here the effects of feeding a (at a logarithmic order, by the still 
large value of its counterpart at the inducing site) and of reducing a by the 
(comparatively small) gel-particle are competitive. 

Let us introduce the generating functions 

m 

Then the evolution equations for Ctmj i™ply the following equation for u, 
which should be valid until the first jump happens (after that u itself will 
be random): 

d d 

—u{x,t,j) = x{u{x,t,j) - Cl^ j - u{l,t,j))—u{x,t,j) 

+ K{u{x,t,l - j) -U{x,t,j)). 

Observe that the first, nonlinear term differs from the one-site situation by 
the appearance of the total mass Ct,oo,j + ^(l>*)i) at the given site instead 
of simply 1 in the one-site expression x{u — l)-^u. The new linear term is 
due to the motion of particles. Of course this equation is not autonomous 
and has to be supplemented by 

which is just the rewritten evolution equation of Ci^oo.j fro™ above. Obvi- 
ously, this system of PDEs looks much less promising for having a (semi-) 
explicit solution, as was the case in the one-site model. 

We emphasize again that all these considerations are heuristical though 
likely to be valid, until a rigorous proof of the convergence to the limiting 
process is available. 

We defined gelation for our model in a weak sense meaning that a tends 
to infinity, which is strong evidence (but of course no rigorous proof) for 
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the simultaneous appearance of gel-particles at the first gelation instance if 
K > 0. In fact, the possible later gelation events cannot be induced gelations, 
since the presence of a gel-particle at the other site keeps the second moment 
there at finite values. So the first gelation is particularly interesting (and the 
proof of the rigorous limit theorem can be expected to be difficult mainly 
because of this special circumstance). 

Bear in mind that, for N being large but finite, only a tiny portion of 
logarithmic order of the quantity a feeds the induced gelation at the other 
site, or, to say it in different words, there is only a small time window around 
gelation time where many very large particles aggregate at the site where 
the "genuine" gelation just started, so only a small part of them can be 
expected to migrate to the other site to induce gelation there. In fact, it is 
quite likely that just as in the mean- field model the largest particles (except 
for the gel-particle) are of size < A^^/^ and these "mesoscopic" particles are 
present only for a short time before and after gelation. They carry the growth 
of the a. The respective part of them which succeeds in jumping to the 
other site during their short time of existence is likely to be of logarithmic 
order in view of the calculations made in the proof of Theorem 4. This 
leads to the situation that, in contrast to the inducing site's growth of the 
expected particle mass per monomer to "infinity" (in fact, some power A^", 
a < |), the growth of the expected particle mass at the induced site is carried 
only by the few very large particles which immigrated and will probably be 
of logarithmic order in A^. Therefore three effects should be likely to be 
observed even for very large of the order of Loschmidt's constant or even 
of the order 10^°°: 

(1) the induced gelation happens with a significant delay of the order 
(log AT)-!; 

(2) since it is carried by only few large "foreign" particles, in contrast 
to the inducing gelation the growth of a in the (delay) pre-gelation interval 
can be expected to show visible erratic deviations (caused by jumps from/to 
the other site and coalescences between the large particles) from the law of 
large numbers; and 

(3) since the initial growth speed of the gel-particle is governed by the 
value of fj at gelation time, which is maximally of logarithmic order in A^ at 
the induced site, the second gel-particle's growth will be slow [we referred in 
Section 3 to the fact that in the mean-field case and hence probably also at 
the inducing site the (unrescaled) primary slope of the gel-particle's mass is 
of order 2A^]. 

Though at the time being no quantum computer was available to us to 
study simulations with lO^*"^ particles (or only 10^^), the simulations with 
A'^ = 10^ rather clearly show the predicted peculiarities of the induced gela- 
tion; see Figure 4 below. In this simulation, the behavior at site is clearly 
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governed by the law of large numbers, except for the random jump of the 
giant particle from 1 to at time 5.86. At site 1 the boost of the second 
moment induced by site around the (inducing) gelation time is visible, but 
not too strong. The delay between inducing and induced gelation is about 
one time unit, and during this delay time the second moment at site 1 is 
mainly carried by just a few particles, which fact causes strong fluctuations 
visible as some large jumps, respectively a "Brownian-like" behavior around 
the (induced) gelation. As criterion for gelation a threshold of size 2 • 10*^ for 
the largest particle was chosen, that is, 2A^^/^. 



5. Remarks concerning other gelling kernels. The situation for general 
gelling kernels (instead of mn) is more complicated, even in the mean-field 
case. Existence of gelling solutions has been conjectured for a long time (see 
the discussion in [1]). Rigorous results were given in [6, 10]. The asymptotic 
behavior of the stochastic model was studied, for example, in [3, 7, 10, 14, 15]. 
Convergence has been proved in the pre-gelation interval. After gelation, 
the Smoluchowski equation might have to be modified (dependent on the 
kernel). Convergence is established under the assumption that uniqueness 
holds. The asymptotic size of the largest particle is still an open problem. 
It has been studied in [1, 16] (see also the corresponding sections in [3, 19]). 

Here we add some heuristic comments concerning the particle behavior 
around gelation, explaining the peculiarity of the mn-kernel in comparison 
to several other gelling kernels. 

Kernel (mn)" for 1/2 < a < 1. If there was a giant particle of size 7A^, 
it would absorb the particles of fixed size m at a rate 



N~\jmNY 



7 m 



that is, at an asymptotically vanishing rate. So it looks a bit surprising that 
it occurs at all. The explanation should be some kind of food chain. Just as 
in the a = 1 case, in the pre-gelation phase the coagulation process shifts the 
total mass progressively to particles of larger size (individual mass), letting 
at the gelation point the mass flow "into infinity" reach an asymptotically 
positive speed. That means that the speed at which mass is lost from the 
set of all particles of size m < B no longer goes to zero as B tends to 
infinity (at a suitable slower speed than N). Hence a positive part of the 
mass accumulates at particles of a size which grows to infinity as N goes 
to infinity. An indicator of that situation is that the mean particle size per 
monomer ("second moment" in the usual terminology) a grows to infinity. 
It is not a priori obvious what should be the typical size of the particles 
carrying that lost mass. 
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Fig. 4. Simulation with one billion particles, A = 0,k = 0.29. 
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Anyway, it might be conjectured that just as in the a = 1 case the "lost 
mass" changes the regime abruptly by bringing the second moment back to 
finite values and stopping the "active" transfer (from the point of view of 
bounded-size particles) of mass to infinity. If this were true, the lost mass 
would actively absorb mass from the finite-size particles at positive speed, 
preferring large particles (that is the reason why the second moment drops 
to finite values for a = 1; it should be noted that here "finite" always means: 
remaining bounded for N tending to infinity, and does not refer to a limit 
model). As was mentioned above, the absorbing power of the giant particle 
is not enough to do this. For a < 1 it must be taken into account that, while 
the rate at which a single particle of size absorbs m-particles is given by 
l^i^a ; for the same mass being distributed on 'yN'^~^ particles of size 
the absorbing force is given by 

Ar-i^iVi-/3(^iV/3)« = ^^«Ar/?("-i). 

The exponent is still negative, indicating that, though splitting the mass into 
many smaller particles increases the absorbing force, the totality of "infinite- 
size" particles is unable to actively withdraw mass from the bounded region. 

Consequently, gelation will go on in a way very different from the critical 
case a = 1: The bounded-size particles do actively transfer mass to infinity 
even after the gelation point, and this mass is transported through a kind of 
food chain at positive speed to infinity, reaching the giant particle only at 
the end of this chain. Observe that the mass transfer from particles of size 
directly to the giant size particle is given by a rate 

A~^(A^A)" = A^"(i+/3)-i 

which tends to zero unless (3 > — — I. So the food chain should be located 
between exponents and ^ — 1- 

It should be expected that, in contrast to the critical case a = 1, the 
second moment a keeps being infinite after the gelation point. 

As for a corresponding spatial model (with mass- independent jump rate), 
the behavior should be significantly different compared to the critical mn 
case: The portion of mass concentrated on finite-size particles might behave 
deterministically after gelation, too, because of a law of large numbers. This 
should be a kind of deterministic spatial Smoluchowski equation. Since the 
giant particle(s) at the spatial sites is/are unable to directly consume mass 
from below the threshold - — 1, the random jumps performed by them (or 
it) should be without effect on the rest of the mass. 

Kernel m^n + mn" for a <1. The behavior will be similar to the mn- 
case, since the giant particle is able to withdraw mass from the bounded 
region at positive speed in this case. 
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Kernel m" + for a > 1. This kernel can be expected to be instanta- 
neously and completely gelling [1, 11]. In fact: For large N after an arbitrary 
short time an arbitrary large particle of size A^' is present with arbitrary large 
probability (of course at the beginning 1 <C A^' <C A'^). This particle by the 
law of large numbers grows like the solution of the ODE x = x", starting at 
initial value N' . So it takes an arbitrary short time to consume the rest of 
the population. Hence no interesting spatial model should be expected. 

Kernel m^n^ + m^n'^ for a + (3> 1. If one of the parameters is greater 
than 1, the argument of the case jtt," + n", a > 1 applies. Otherwise, if a, /3 < 
1, it seems that anything said about the (mn)" case could be valid here, too. 

6. Proofs. We start by quoting a result on differential inequalities in 
which will be used in the following considerations at several places. Let for 
X, y G R*^ the notation x <y mean that Xj < yi for i = 1, 2, . . . , d. A function 
called quasi-monotone increasing (with respect to this natural 
semi-ordering) if for each z G M*^, < 2 from x < y and {z,x) = {z,y) it 
follows that {z,g{x)) < {z,g{y)). In fact, we shall employ this definition only 
for the cases d= 1 (where any function is quasi-monotone increasing) and 
d = 2, where it simply means that gi{x) is a monotone increasing function 
of X3_j when Xj is fixed, for z = 1,2. Then we have (cf. [9] or [17]): 

Theorem 5. Let f:[a,b] x M*^ ^ M'' be a function such that fit,-) is 
quasi-monotone increasing and locally Lipschitz continuous for each t € 
[a,b], and let u,v:[a,b] -^M."^ be dijjerentiable. Then 

u{t) - f{t,u{t)) < i{t) - f{t,v{t)), t G [0,6] 

and 

u{a) < v{a) 

imply 

u{t)<v{t), t£[a,b]. 

6.1. Proof of Theorem 3. First observe that the theorem about existence 
and uniqueness of ODE ensures that there is a unique maximal local smooth 
solution to this evolution equation up to some < Tmax < +00, and if T^ax < 
+00, then necessarily the vector {x,y) becomes unbounded near Tmax- Con- 
sider S := X -\- y. It fulfills the equation S = x^ -\- y"^ > ^(x -|- y)^ = ^S^. 
The maximal solution of the corresponding ODE is (1 — ^t)^^, <t < 2. 
Hence we see (using Theorem 5 in the case d = 1) that Tmax < 2 is fi- 
nite. Let tgei := Tmax- We know that at least one of the x, y will be un- 
bounded near t^^i- Assume it is the function x. Since x > x^ — kx > — kx and 
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y > y'^ — Ky > —Ky , the functions are nonnegative for t>0 (again by Theo- 
rem 5). Furthermore, since x is unbounded near igei, this estimate combined 
with Theorem 5 shows that it is in fact monotonically increasing near tgei 
and hence tends to +oo. Assume y is bounded near tgd. Then there is some 5 
such that y"^ + k{x — y)>^KxfoTt> tgd — 6 implying y > ^kx for t > tgd — S. 

Now represent x as g{tgci — in a left neighborhood of tgei with g being 
smooth in a right neighborhood (0, 6) of where we may use the same 5 as 
above. We get 

x = x"^ + K{y - x) = g{tgci - t)~^g(tgcl - t) = x^gitgd - t), 

so that 

g(t) = [{x{tgcl - t)f + K(y(tgel - t) - x{tgci - t))]/ (x(tgcl - t)f 1 

as t — > by assumption. Hence g{t) = c{t)t with c being bounded from above 
and below by a positive constant C (resp. its inverse C~^) in a right neigh- 
borhood of 0. We derive 

y{tge\ -t)= y(tgci - 5) + / y{u) du > y(tgci -S) + - / x{u) du 

1 rtgci-t 

= y{tgc\ - 5) + -K c{tgci-u) (tgci-u) du 

> -kC I (tgei - du = -KC(ln(5) - Int). 

The last expression is unbounded as f — > in contradiction to the assump- 
tion. Hence we know that y is unbounded near tgd, too. As above we may 
conclude that y — > +00 as t ^ tgei monotonically. We have shown the asser- 
tion. 



6.2. Proof of Theorem 4. Consider D = x — y, S = x + y and assume 
3^0 > yo- These functions fulfill the quadratic ODE 

5 = i52 + iZ)^ D = SD-2kD 

with initial condition = 1, < Dq < Sq. Since S is obviously monotoni- 
cally increasing we get D = SD — 2kD > (1 — 2k) D from which we infer by 
Theorem 5 that Z? > is always fulfilled. Hence the site with originally more 
monomers has always a larger expected particle size. And, since S goes to 
infinity monotonically near tg^, we conclude that D is positive near tg^i, so 
D increases monotonically near tgei, too. Assume D is bounded 
above, represent S as /(tgci — 0"^ ^ 1®^^ neighborhood of tgd which yields 

S=^S^ + = /(igel - t)-^f{tgd -t) = S^fitgd - t) 
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and we may infer again that / would tend to ^ near and as above we 
would obtain a contradiction, since in this case the integral 

S{u)D{u)du>D{t^^i-5) / ° _ S{u)du 



diverges near tg^. So in fact D increases monotonically to infinity near tgd. 
Let Q := f < 1. Then it is easy to see that Q fulfills 

(14) Q = \D{1-Q'^)-2kQ. 

If Q was bounded away from 1 near tgei we might conclude from this differ- 
ential equation that it is strictly increasing near tgd, hence having a limit 
Qo < 1- With 

S = \S^ + = /(tgel - t)-V(igel -t) = 5V(tgel " t) 

we would conclude that / would tend to ^{1 + Qq) near and hence 
J^^''l_gS{u) du as well as JIi^°I_^ D{u) du would diverge near tgd- This would 
imply the divergence of Q near tgd in view of (14), a contradiction. We have 
shown that 1 is an accumulation point of Q near tgc\. If it was not the limit, 
Q would necessarily assume arbitrary large negative values in contradiction 
to the boundedness from below by —2kQ. We have shown that Q ^ 1 as 
t tgci- This means that y/x —i-O as t —>■ tgd, provided xq> i/q. Repeating 
the calculations some lines above, this time with the confirmed value Qq = 1, 
we get that 5 = /(tgd — where / — > 1 as t — > 0. So there are smooth (in 
a right neighborhood of — and excluding — 0) functions c, ci,C2 all tending to 
1 as i — > such that 

S = Ci(tgel - t){tgel - t)"\ D = C2(tgel - t)(tgel " 

X = c(tgel-t)(tgcl-t)"\ 

if xq > uo- It remains to calculate the asymptotics of y. 

We consider the quantity q :=y/x < 1 which fulfills q = q'^ {x — k) — qx + k. 
We know that q is positive but tends to zero near tgei, and qx = y tends to 
infinity. Hence for t sufficiently close to tgd we have q < —(1 — 6)qx with 
arbitrary small 6. So 

[logq] < -(1 - s)x < -(1 - ,5)c(tgd - t)(tgci - ty^ < -(1 - 2<5)(tgd - ty^ 

for t sufficiently close to tgd. Consequently we get logq < cq + (1 — 25) x 
log(tgei — t) (for some cq) in a neighborhood of tgd- So q < e'^°(tgd — t)^""^^ 
and hence y < e'^''(tgei — t)"^^ . This proves that y is not only o{x) but even 
o{x^) for each e > 0. 

We return to the original equation y = y"^ + k{x — y). Now it is obvious 
that for ?/o < ^^o the increase of y near tgd is (asymptotically) completely 
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caused by the migration term kx, while the coagulation term can be 
neglected. For any positive 5 we find for t sufficiently close to tgci: (1 — 5)kx < 
y < {1 + 6)nx and the x-asymptotics yields 

-(1 - 2(5)Klog(tgel -t)<ko + y<-{l + 25)Klog(tgei - t) 

in a sufficiently small neighborhood of tgd, where fco is a constant de- 
pending on 6. Therefore we finally deduced the logarithmic asymptotics 
of y = — C3(tgei — t)Klog(tgci — t) with C3 tending to 1 as t ^ 0. 

6.3. Proof of Theorem 2. 

Existence. Fix some T, < T < t^^^. Consider the following system of 
ODEs: 

^t,m,j ~ 2 5Z ~ ^)^t,n,j^t,ra—n,j ~ ^'^t,m,jCt,*,j 

0<n<m 

with initial condition 

•^0,1,0 — Co,*,o ~ ]^ _|_ _>^' '^0,1,1 — Co,*,i ~ -j^ _(_ ^' m J = 0, m>l. 

It has obviously a unique solution for all t since, given i? G N, the equations 
for m < B, together with the equations for (^ ^^j form a finite autonomous 
system. First observe that all quantities Cf m,j}Ct,*,j nonnegative. In fact, 
this is easily seen for (^^^ j which can be given in explicit terms: 

To see it for the c^ ^^j, let us first modify the initial condition to Cq ^j — 
c, m > 1, where c > is an arbitrary fixed value. It is enough to prove the 
nonnegativity of the corresponding solutions for all times t > in these 
cases in view of the theorem about the continuous dependence of solutions 
of ODEs on the initial data. Assume the opposite. Then there would be 
some c such that at least one of the c^^j reaches for the first time at 
some time < to- We get the differential inequalities c^ mj ^ {—m — K)c^ .^ j 
being valid for < t < to taking into account (16). But now we immediately 
infer from Theorem 5 that c^^ ^ j > in contradiction to our assumption. 

Hence indeed all c^^^ are nonnegative for all times. 
Next, consider the quantities 

m<B 
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We have 

^t,*,j ~ h ^ ~^ ^)^^'^t,k,j'^t,l,j ~ Ct,*,j 5Z ''^^^t,m,j 
k+KB m<B 

_i_ (r^^^ —/■^^^\ 

(17) 

< ^ /C lCt^k,j(^t,l,j - Ct,*,j X! "^t.mj + '^(Cj,*,l-j - Ct,*j) 
k,l<B m<B 

which imphes 

The initial values 5^'^ are zero. So we may conclude, again by Theorem 5, 
that 5^ j is nonpositive yielding 

(18) c5j<Ct^j<i. 

Next, consider the quantity 

B,X 2 A 

m<B 

We have 

k+l<B m<B 
, I B,X B,X^ 

(19) < Y: kHcl,^,cl^^ + (ag^)2 - d^.^ ^'<m, 

k,l<B m<B 
, / B,X B,X^ 



/ B,X\2 , /^B,X y-X \ 3 A , / B,X B,X\ 

m<B 

Hence, by the above estimate Q^i^j < Ct,*j arrive at 

■ B,X ^ , B,X^2 I / B,X B,X^ 

The function defined by fj{x) := x'j + K{x^-j — Xj) is quasi-monotone in- 
creasing. Consequently by Theorem 5, ((j^^j''*)j=o,i is dominated in [0,tgpi) 
by the solution of the ODE (5)-(6), for all B. So (cr^j)j=o,i is bounded 
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from above by this solution, too. Hence a^j is integrable in any compact 

subinterval of [0,tggi). 

Next, consider again the relation (17) which implies 

I k+KB 

m<B ) 

The quantity | Z]fc+«<_B(^ + O^^^^fc j^^/ j converges monotonically (as B —)■ 
oo) toward 



_ .A 
kl m 



for each s € [0,tgpi), in the same way J2m<B''T^'^(^s,ni,j converges toward a^ j, 
and J toward 

m 

[by (18)]. So from Lebesgue's theorem about monotone convergence we in- 
fer, taking into account the uniform boundedness of all these quantities in 

[o,r], r<t^,i, 

(20) C;' = C,5 + /J{<,(C4 - ct.) + '^(Ci-. - Cj)} ds. 

Let := - < [by (18)] and = Id^^s' + Q'\ = "Co' " C 
and we get [bearing in mind that for each s we have Cs'* o + Ci"* i = 1 i^i view 
of (15)] 

(21) = / V^oCd' - <iQ'} ds<C f d^ ds 

Jo Jo 

for some constant C depending on T, since a^j is uniformly bounded for 
< s < T < tgp[. Since d^ is bounded by 2, we get from (21) for s < that 
dg is bounded by 1 and by induction that it must be zero in [0, ^] , which, 
arguing by repetition, leads to the conclusion that d^ is zero everywhere in 
[0,tggj). Hence we have = Cs,*,j everywhere in [0,iggj). Now, substituting 
~ ^m''^'^s,m,j ^s,*,j ™ (15) that the unique solution of (15) 

also satisfies (2) in [0,tggj). We proved that this solution is positive and 
strong in the sense of integrability of a^j in any compact interval contained 
in [0,di). 
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Uniqueness. Let a strong positive solution (c^^mj)meNj=o,i,tG[o,t )' ^ ^ *geb 
of (2) be given. Using again the notation (^^j := J2m<B '"^'^t,m,j ^^^^ C^^'j' := 
Yl,m'^^'^t,m from (2) we infer in the same way as in (17) 



^B,\ ^ ,^B,\ _ .COM 2 A , K(t^'^ - t^'^'' 

m<B 



which in view of Cfjj - ^t^*J y^^^"^^ 4^* J - ~ '^Mj) ^^'^ conse- 

' B X ' B X B X 

quently Ct * o + Ct * i ^ 0. Thus the sum of the two positive quantities Ct * j^j — 
0, 1 is always less than or equal to its initial value 1 (for B > 1). Letting B 
tend to infinity we arrive at the conclusion that both Ct^'o'Ct^'i 
greater than 1 for <t <t. Now we have 

•'^ k+KB 

rt rt 

- ^'^^s,m,jds+ fi{Cs,*A-j-Cs,*))ds. 

•^0 m<B •'^ 

With B ^ oo the first integral on the right-hand side tends toward 

f \Y.ik + l)klclk^jcli^^ds= faljQjds< faljds<+oo 

hi 

for each < t < t by the assumption that the solution is strong. The same 
argument yields the finiteness of the second integral, which tends to the 
same limit Jq a^jC^^j ds . So we arrive at the system of two equations 

— ^0,*,j + '^\^s,*,l-j ~ ^s,*,j) "-^^ J — "Ji J-) 

from which we may easily conclude that the bounded quantities j 
continuous and even continuously differentiable and fulfill the ODE 



^oo,A *oo,A I ^oo,A 



with initial condition 



^oo,A aCxDjA 

'50,*,0 ~ 1 I \ ' ^0,*,1 



1 + A' "^•*'^ 1 + A 

Therefore our solution also fulfills (15) which has a unique solution for all t. 
So the maximal strong positive solution must be unique. 
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Maximality. We still have to show that the maximality implies i = tg^^j . 
Knowing already that the solution to (15) in [0,tgj,j) is positive and strong 
[in view of the fact derived above that a^j is bounded from above by the 
continuous solution of (4)], we just need to show that t cannot be larger 
than tgp[. For this sake we first prove that the second moments ((T^j)j=o,i, as 

announced above, fulfill (4) for < t < tg^j. Denote by pfj^,p^ j the quantities 
J2m.<B ''^^'^t,m,j^J2m "^^c^.mji respectively, given the solution to (15). We get, 
in an analogous way as in (19) for 0<t <r<tL, 

k+KB m<B 

, / B,X BM 
+ t^{Pt,{-j - Pt,j ) 

k,l<B k,l<B m<B 

, / B,X B,X^ 
+ l^iPtA^i - Pt. 



3 ^t,J 

BM 



o B,X B,X I f/-B,X aA \ 4 A , / B,A B 



m<B 



^ o B ^X I / B^X B^X\ 

< Sp^ 'j C + K^KPtx^^ - PtJ ) 



where C is a constant depending on T, where we made use of the bounded- 
ness of a^ j in compact subintervals of [0,tggi). Prom this, in the same way 

as for the second moment, we derive that all p^j^ and hence also p^j are 
bounded by the solution of the coupled system of ODEs 

X = 3Cx + K{y — x), y = 3Cy + K{x — y) 

with xq = j^,yo = in [0,7"]. Now we may rewrite the first equation in 
(19): 

E {k + lfklc^s,kAl,jds 



•"^ k+l<B 
/"Va B,X j , /"S B,X B,X^, 



With S — > oo the first integral on the right-hand side tends monotonically 
toward 



Jo ^ Jo 



Pi, + {^sf?)ds 
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being finite for t < tggj, and the second integral converges monotonically 
toward 

wliich is finite for these t, too. So we may pass to the limit on both sides 
and get 

Thus the a^j are continuous and even continuously differentiable in [0,iggj) 
and fulfill (4). Now from Theorem 4 we easily infer that the condition that 
a^j be integrable for each compact subinterval of [0,t) (strong solution) 
implies i < t^^^ . 

6.4. Proof of Theorem 1. We are now prepared to prove Theorem 1. At 
this point we take the point of view that the configuration at time t of our 
particle system is given by a finite counting measure on N x {0, 1}, 

(22) Xt = ^ ^t,m,jS{m,j) ■ 

Let us, just for the moment, renounce the assumption that the initial con- 
figuration be monomeric; that is, Co,m,j niay be positive for m> 1. 

In order to prove our limit theorem we consider a coupled system of 
truncated versions of the process {Xf} indexed by i?EN = NU{oo} where 
the dynamics of {X^} are given as follows: 

1. First define 

Xq ■■= lm<B^(mA)XQ{d{m,i)) + I lm>BmSfB,o)Xo{d{m, i)) , 

JNx{0,l} JNx{0,l} ~ 

that is, the particles of sizes less than B remain unchanged, while at each 
site the masses of all particles larger than B are combined and this mass is 
assigned to (.8,0). Xq can be represented as a finite sum J2i£j ^(miji)^ where 
J is a finite index set, and each Xq can be written as 

(23) J2 \mM+X^i{{B,0)})5^B,o), 

where J{B) is the subset of J corresponding to particles of size smaller than 
B. Obviously X^ = Xo. 

2. It is easily seen that, by simply cutting each Xf in an analogous way, 
the corresponding quantities will not form a Markov process except for k = 
(which corresponds to the one-site case). We intend to let Xq evolve 
Markovian, thereby keeping it in close distance to Xt for B large, that is. 
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to introduce a coupled evolution of {Xt,Xf) where both marginal processes 
are Markovian and their distance can be controlled to be small for large 
B. To perform this, first introduce a new discrete parameter assigned to 
the particles taking two possible values v and h ("visible," "hidden"). That 
means, instead of N x {0, 1} we choose N x {0, 1} x {v, h} as phase space. 
All particles present in Xq with masses below B are marked visible; the rest 
(with masses > B) is marked hidden. This gives X^ . Given a path {Xt)t>Q 
we construct {X^)t>Q as follows, not changing the projection to N x {0,1} 
which is {Xt)t>o-- 

(a) the mark "hidden" is dominant; that is, the combination of two par- 
ticles occurring in {Xt)t>o produces a hidden particle iff at least one of them 
was hidden, and jumping to the other site does not change the status hidden 
or visible, 

(b) if a particle of mass at least B is generated by the combination of two 
particles, it is marked hidden, and 

(c) additionally, particles become "hidden" by a newly introduced inter- 
action between particles of different sites: given X^ , a "visible" particle of 
mass m< B at site j becomes "hidden" at rate 



So the interaction concerns only the newly introduced mark, and has the 
effect that any pair of particles with masses m, k at different sites assigns 
the mark "hidden" to the first particle at rate mk, assuming the second 
particle is already hidden. This represents a Markovian refinement of the 
original process. 

Now we define Xj^ by 



This means that only visible particles are present in X^ at locations < 
whereas all hidden particles (in particular those with masses at least B) 
are collected to form the mass at B at site 0. Observe that the total rate at 
which a visible particle of mass m becomes "hidden" is given by the product 
of m with the total mass of all other hidden particles. 

3. So, if we do not intend to describe the joint dynamics of {Xt,XP) but 
only that of Xf we may do it as follows, referring to the representation (23): 
Consider a family {Ti^ii)i^i'^j(B),i<i' ,ji=ji, of independent exponentially dis- 
tributed random variables with parameter mimi', a second (independent of 
the first) family (5';'^)/ej(s) of independent exponentially distributed random 
variables with parameter kN, and a third independent family {UP)i,zj(^B) 
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of independent exponentially distributed random variables with parameter 
miX^{{{B,0)}). These families are finite in view of the finiteness of J. Let 

T = rn\n{Ti^,,SP,Uh. 
Define = X^ for t < T and 

mi + mil < B, 

mi,mi> < B, 
m +mii > B, 

^0 - \mi,ji) + ^{mi,i-ji), if T = SP,mi < B, 

^0 -^imi,ji)+m\B,0), ifT = Ui^, 

^ Xq , otherwise. 

Starting this construction anew with X^ instead of Xq we finally define X^ 
for arbitrary t>0, thus arriving at a Markov jump process {Xj^}. Observe 

that {Xn = {X'^^'^'^+'j = {Xt}, where u{Xo) := Jj,^^o,i}^Md{m,i)) is 
the total mass. 

The evolution of a single X^ may be described as follows: 

• Particles of sizes m,n at the same site coalesce at rate mn, assuming 
m,n < B. If m + n < i?, a new particle of size m + n is generated, whereas 
for m + n > B the mass m + n is added to the mass already aggregated 
at {B,0). 

• Particles of size less than B change their location at rate kN. 

• Finally, particles of size m < B are deleted at a rate given by the product 
of their size and the mass already aggregated at {B,0) and their mass is 
added to the aggregated mass. 

The evolution conserves the total mass 

u{Xo)=UB{Xi) := / mX^{d{m,i))+X^{{{B,0)}), 
JNx{0,l},m<B 

being independent of B. 

Let, for B Gf^, Mb denote the space of finite signed measures on 

(Nn[i,i?])x{o,i} 

and define -.Mb Mb by 

= \ I [f{x + y,j)lx+y<B 

J(Nx{0,l})2 

+ f{B, 0)(x + y)lx+y>B - fix,j) - f{y,j)] 
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X lx,s/<Bj=fc • xy{n X fi){d[x,j], [y, k]) 



+ / M{(^,0)})[/(i?,0)x-/(x,i)]l,.<BX/x(d[x,j]) 
■/Nx{0,l} 

+ K [f{x, 1 - j) - f{x,j)]lx^BtJ'{d[x,j]), 

JNx{0,1} 

for any bounded function / on N x {0, 1}. Fix a finite measure [Iq on x 
{0, 1} with 



v{^q)= x/io(d[x,j]) < +00 

lx{0,l} 



and define 



^0 := / '^m<B5{m.i)m{d{m,i)) 



+ / lm>Bm8lBfi)^^G{d{m,i)). 

JNx{0,l} 



It is easy to see that for each i3 G N the equation 

(24) = ^^+ fL^{f^f)ds 

Jo 



has a unique solution which is a continuous map fi^y.t —>■ Mb- In fact, 
writing 

i<B,je{0,l} 

results in a finite approximation of the two-site Smoluchowski system: 

k+m=i k<B 



(25) + K(cfi,i-j - cf^i,j) for i < B, 

^t,Bfi = \ X! ''^)^''^'^t,k,j'^t,m,j + C^B,0 X! ^^'^^fcj' 

k+m>B k<B 
k,m<B 

This system of i? first-order ODEs obviously has a unique global solution. 
We may apply Theorem 5 in the same way as was done in the first part of 
the proof of Theorem 2 to see that the solution is nonnegative. Moreover, it 
is not difficult to check that (25) conserves the total mass 

B\ . f B/ 



^im) ^ ^B{^^") ■■= / mf,^{d{m, i)) + fi^{{{B, 0)}) 

JNx{0,l},m<B 

being independent of B. 
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Proposition 6. The solution of the approximating system (25) with 
initial condition 

^0,1,0 = _|_ ^1 = _|_ '^o,m,j =0) m > 1 

fulfills 

B-l 

E ^I^mj - cmjI ^ c^B.o < ^"'c^o, < t < r < 

i=l 

where Cq is a constant depending on T and A. 

Proof. Let us recall the two-site modified Smoluchowski equation (omit- 
ting the upper index A for simplicity) 

(26) . 

which we considered for the initial condition 

C0,1,0 = Co,*,0 = 1 , X 1 Co,l,l = Co,*,l = -.,, > CO,mj = 0, TTi > 1. 

1+A 1+A 

The solution obeys Ct,*,j = J2k ^ct,k,j up to tgei as was shown in the proof of 
Theorem 2. Moreover, using the quantity 



k>B k<B 

we may infer for t < tgd from (26) 

k+m=i k<B 



(27) 



Define 



- + KCt,i,i-j lor i<B, t< tgel, 

k+m>B k<B 
k,m<B 



'■— '^t,i,j C^jj) i<B. 
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Subtracting (25) from (27) we get 

k+m=i k+m=i 
k<B k<B 

k+m=i 
k<B 

The mass conservation law for (25) and (27) implies 

^^B.o = 1 - E E ^'4-,^ = + dfl, + E E - ^^m)- 

1=0 k<B 1=0 k<B 

Thus we obtain 

k+m=i 

— i E ^iPu,j'^t,k,j + Ct^i,jPt^k,j) 



k<B 



~ ( + E ^^(^^k,j + Pt^k,l~j) I 

V \ fc<B / 



fe+m=i 

- ^ E ^f^t!i,jCt,k,j 

k<B 

Finally we arrive at the following finite system of ODEs: 



k+m=i 



\ k<B J \ k<E 



B 

3 

k<B 
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with initial condition /3(fj j = 0, which has a unique solution. Again we may 
apply Theorem 5 in the same way as was done in the first part of the proof 
of Theorem 2 to see that the solution is nonnegative, taking into account 
the boundedness of the total mass of the solution to (27) by 1 and the fact 
that this solution was already proved to be nonnegative. Hence we have the 
relation 

(28) cf^^J<ct,i,,, i<B,t>0. 

Thus 

B-l B-l 

\^t,i,j — '^t,i,j \ = X/ ~ ^t,i,j) 

i=l i=l 

and since both systems conserve the total mass we get 

1=1 

So, in order to estimate the total deviation of the solution of (25) with 
respect to that of (26) for masses less than B it is enough to estimate q. 
We get from (25) and (28) for any t < T < tg^i 

k+m>B k<B 
k,m<B 

^ I + '^)^^Ct,k,jCt,m,j + C^Bfi ^'^'^tAj 

k+m>B k<B 
k,m<B 

^\ Y + ^)^^<^t,k,jCt,m,j + Ct^B,oic^tfi + Crt,l) 
k+m>B 
k,m<B 

^ I Y + ^)^^'^t,k,jCt,m,j + Ct^BfiC, 
k+m>B 
k,m<B 

where C is a constant depending on T as it was shown in the proof of 
Theorem 2 that the second moments at,j = cr^j are continuous functions 
in each compact subinterval of [0,tgei). Hence in [0,T] we have (using the 
notation [x] for the integer part of a real number x, and Ci,*,o + Ci,*,i = 1 for 
the fourth line) 

cfBfi ^ Y^ k'^'^Ct^k,jCt,m,j + (^f^Bfi^ 
k+m>B 
k,m<B 
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< k'^rnct,k,jCt,m,j + cf^BfiC 

k+m>B 

k>[B/2] rn>[B/2] 
m k 

< k^ct,k,j + (0-4,0 + o-t,i) Y '^'^t,m,j + cIb,oC 

k>[B/2] m>[B/2] 

< Y f^^Ct,k,j + C Y "T-Ct.mj + Ci^B.oC 

k>[B/2] m>[B/2] 

<[B/2]-4 k''ct,k,j+C Y ^^^t,m,j]+cf^B,oC 

\k>[B/2] m>[B/2] I 

< [B/2]-^ k^ct,k,j + C{at,o + at,i)j + cf^BfiC 

< [i?/2]-i (^Y f^'^tAj + C'^ + cf,B,oC 

<[B/2]-HC' + C') + cf^^^,C 

as, again in the proof of Theorem 2, it was shown that also the third moment 
is bounded uniformly in compact subintervals of [0,tgei). Hence we see that 
the positive quantity cf^ q in [0, T] is bounded (in view of Theorem 5) by 
the solution of the differential equation 

x = [B/2]~'^{C' + C^)+Cx 

with initial condition x(0) = 0, which is given by 

x{t) = [B/2]^^C^^{C' + C^){e^' - 1). 

Therefore, there is a constant Co depending on T such that cf^ q < B^^Cq 
in [0,r]. This proves the proposition. □ 

Consider the process (22) and the corresponding coupled family of ap- 
proximations constructed above. Set 

Xt = N''^Xf/jM, = N~^X^i^. 

Let d be a metric on the set of all finite measures on N x {0, 1} which 
generates the weak topology, that is, (/, ^n) {f, lA if and only if d(/x„, /x) — > 
for each bounded function on N x {0, 1}. We may choose d in a way that 
din,fJ'') < WfJ' — /i'll for each pair of measures, where || • || denotes the total 
variation norm. 
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Proposition 7. Assume 
(29) No/N^j^, ^i/^^Y^ asiV^oo. 

Then, for all t>0, 

sup(i(Xf,//f)^0 

s<t 

in probability as N ^ oo. 

Proof. The sequence of laws corresponding to {Xf)t>o is tight, and any 
weak hmit point (with respect to N) is concentrated on the set of solutions 
of the equation (24) (cf., e.g., [4], Theorem 2.3). Since the solution (/if) of 
this equation is unique, we get convergence in distribution, which implies 
the statement of the proposition (cf., e.g., [18], Section 3). □ 

The following theorem and its proof adapt Theorem 4.4 of [14] to our 
situation. It is a stronger version of Theorem 1. 

Theorem 8. Assume (29) and consider the solution (/Lft)o<t<t^ '^f 
Smoluchowski equation up to tg^j. Then for all t < t^^j 

supd{xXs{d[x,j]),xfis{d[x,j])) 

in probability, as N ^ oo. 

Proof. For simplicity, we omit the upper index A and write dx instead 
of (i[a;,j]. Fix 5 and t < igei. By Proposition 6 we find some B such that 
/if ({(i?, 0)}) < S/A. Since the measures are supported on a compact set. 
Proposition 7 implies 

supd{xX^{dxn[l,B)),xfj.f{dxn[l,B)))^0 

s<t 

in probability as — > oo. Using that the derivative of /zf ({(S, 0)}) is non- 
negative, we have by Proposition 6 for s <t 

\\xfis{dx) - x^if{dx n [1, B))\\ < 2fif{{{B, 0)}) < 2/if ({(S, 0)}) < S/2. 

On the other hand, 

\\xXs{dx) - xX^{dxn[l,B))\\ <2X^{{{B,0)}) 

by the fact that the evolutions of particles smaller than B coincide for {Xf} 
and {X^} except for those which are hidden and aggregated at (-6,0) in 
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{X^}, and by conservation of the total mass. We may proceed as follows, 
bearing in mind that the aggregated mass is increasing in time: 

\\xXs{dx) - xX^ {dx n[l, B))\\ 
<2Xf{{{B,0)}) 

<2/.f({(i?,0)}) + 2|Xf({(i?,0)})-/.f({(i?,0)})| 
<S/2 + 2\\Xf - 

So, using the property of d to be dominated by the variation distance, we 
obtain 

d{xXs{dx) , xfj,s{dx)) 

< \\xXs{dx) -xX^{dxn[l,B))\\ 

+ d{xX^{dxr\[l,B)),xfif{dxr\[l,B))) 

+ \\xi2s{dx) -xfifidxn [l,B))\\ 
<5 + 2\\X^ - H^\\+d{xX^{dxn[l,B)),x^if{dxn[l,B))) 

which gives 

pi sup d{xXs{dx),x^s{dx)) > 6] -^0 

\s<t / 

by Proposition 7, bearing in mind that all metrics considered here are equiv- 
alent on the finite-dimensional vector space corresponding to measures on 
{l,2,...,S}x{0,l}. □ 
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